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Abstract 



The two matrix model is considered with measure given by the exponential of a sum of 
polynomials in two different variables, ft is shown how to derive a sequence of pairs of "dual" 
finite size systems of ODEs for the corresponding biorthonormal polynomials. An inverse 
theorem is proved showing how to reconstruct such measures from pairs of semi-infinite finite 
band matrices defining the recursion relations and satisfying the string equation. A proof is 
given in the N — > 00 limit that the dual systems obtained share the same spectral curve. 

1 Introduction 

We consider the two-matrix model [0, 0, 0, 0], which involves an ensemble 
consisting of pairs of N x N hermitian matrices Mi and M 2 , with a U(N) 
invariant probability measure of the form: 

—dfj,(M 1 ,M 2 ) := — exp/ftr {-V X {M X ) - V 2 (M 2 ) + M 1 M 2 )dMidM 2 ■ (1-1) 

TN TN 

Here dM\dM 2 is the standard Lebesgue measure for pairs of Hermitian matrices 
and Vi and V 2 are chosen to be polynomials of degrees di + 1, d 2 + 1 respectively, 
and are referred to as the potentials. The overall positive scaling factor K in the 
exponential is taken as having order N when considering the large N limit. We 
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2 bertola@crm.umontreal.ca 
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also assume that both potentials are real and bounded from below (for reasons 
of convergence) . 

The normalization factor (partition function) 



r N = dfx (1-2) 

J Mi JM 2 

is known to be a KP r-function in each set of deformation parameters (the 
coefficients of the two polynomials Vt, V2), as well as providing solutions to the 
two-Toda equations |n], [I], ||. The key objects of the theory are the correlation 
functions for the eigenvalues of the two matrices. Analogously to the one- 
matrix models, such correlation functions can be recovered by means of certain 
Fredholm integral kernels. We recall here briefly that in one-matrix models with 
measure 

—dfi(M) := — exptr (-V(M))dM (1-3) 
tn t n 

the computation of such a kernel is reduced to the construction of orthonormal 
polynomials P n (x) for the space L 2 (M.,e~ v ^dx) . In terms of these polynomi- 
als, the kernel is given by 

N-l 

K(x,x') = Pn(x)e-^P n (x')e-i v ^ . (i_ 4 ) 

n=0 

In 2-matrix models there are four relevant kernels needed to compute the sta- 
tistical correlations of eigenvalues. For m-matrix models there are m 2 such 
kernels. 

These kernels are expressible in terms of suitably defined sequences of biorthog- 
onal polynomials. By this we mean two sequences of monic polynomials 

n n (x) =£" + ••., a n (y) =y n + ■■■, n = 0,1,... (1-5) 

which are orthogonal with respect to a coupled measure on the product space: 

[dxdy n n (x)a m (y)e- KV ^- KV ^ +Kx y = h n S rnn , (1-6) 

R 

where V\ (x) and V% (y) are the polynomials appearing in the two-matrix model 
measure (1-1). The orthogonality relations determine the two families uniquely, 
if they exist ||. The four relevant kernels are expressed as follows in terms of 
these biorthogonal polynomials 



N Ar_1 1 

Kv2{x,y) = T-*n{x)a n (y)e- KV ^e- KV ^ , (1-7) 



n=0 
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Ku(x,x') = J dy Ki2(x,y)e Kx 'y, (1-8) 
K22(y',y) = [ dx Ki2(x,y)e K *v' , (1-9) 

JR 

n r r n , , 

K21 (y',x')= dx dyKi2 (x, y) e Kxv e Kx » . (1-10) 

All the statistical properties of the spectra of the 2-matrix ensemble may then be 
expressed in terms of these kernels Q and the corresponding Fredholm integral 

N 

operators J&j, i,j = 1,2. For instance the density of eigenvalues of the first 
matrix is: 

N IN 

Pi{x) = — Ku(x,x) , (1-11) 

the correlation function of two eigenvalues of the first matrix is: 

jv , 1 / N N N , N , \ 

Pii(x,x') = j^iKii(x,x)Kii{x',x')-Kii(x,x')Kii(x',x)j , (1-12) 

and the correlation function of two eigenvalues, one of the first matrix and one 
of the second is: 

N 1 (N N N N \ 

Pi2(x,y) = j^iKii(x,x)K22(y,y)-Ki2(x,y)(K2i{y,x)-e Kxv )J . (1-13) 

Any other correlation function of m eigenvalues can similarly be written as a 
determinant involving these four kernels. 

The main objective of this paper is to derive and analyze certain differential 
systems of ODE's satisfied by the quasipolynomials ipnix) := n n (x)e~ Vl< -- x \ 
0n(y) : — o'niy)^^ 2 ^ and their Fourier Laplace transforms. In section [| we 
summarize the principal results for finite TV. The details and proofs may be 
found in (|] and Q. In section || we derive the corresponding results in the 
N — > 00 limit in a simple way. 



The proof of Prop. 2.1 and the non-abelian version of the transvcrsality 
argument in Section [| is based on joint work with J. Hurtubise, details of which 
will appear in Q. 

2 Folding and systems of ODE in duality 

Consider the normalized quasi-polynomials 

1> n (x) - -^=Tr n (x)e- KV ^ , 0„( y ) = _L a n (y)e- KV *M , n = 0,...oo . 

(2-14) 
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Viewing these as components of a pair of semi-infinite column vectors 

* = (ip ,tp 1 ,...ip n ,...) t and $ = (0o,0i, ...<t> n ,---) t , (2-15) 

oo oo 

we obtain a pair of semi-infinite matrices Q and P that implement multiplication 
of by x and derivation —-kS:, respectively. Equivalently, we obtain the 

oo 

transposes Q l and P* by applying — -h^j or multiplication by —y to $. By 
construction, these satisfy the Heisenberg commutation relations (or "string 
equation" ) 



[P,Q] 



1 . 

K 



(2-16) 



Along with these quasipolynomials we need their Fourier-Laplace transforms 
and the corresponding semi-infinite (row)-vectors with components 

dye K ^4, n (y) (2-17) 
&,...,£,,...)• (2-18) 



±Jy) := / dxe K *"M*) > ±A*) 



The multiplicative and derivative recursion relations for these sequences can be 
shown (by integration by parts) to be 



oo oo GfX OO oo 

v*{v) = Mv)Q t ; 4^-*(f) = Mv)P t ■ 

oo oo " Oty oo oo 



(2-19) 
(2-20) 



It also follows |3(] from integration by parts that the two matrices P and Q have 
a finite band structure 

a (0) 7 (0) ..." 

ai(l) ao(l) 7(1) ••• 



Q:= 



P := 



a (d 2 ) 7(^2) 



0^2(^2) 

'■• 

A)(0) /3i(l) 

7(0) /3 (1) ft (2) 

7(1) A)(2) 

7 (2) /3 (3) 



(2-21) 



^(di + l) 



(2-22) 



where 7(71) 7^ for all n E N. This structure essentially follows from the fact 
that the two matrices 

(P-V{(Q)), (Q-VUP)) (2-23) 
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are strictly lower and upper triangular respectively. Indeed, in the basis of 
quasipolynomials it is obvious that 

£ (P - V{(Q)) nm ip m (x) = (~~ ~ V{(x)) Mx) 

m=0 ^ ' 

= cip n -\ (x) + lower degree quasipolynomials. (2-24) 

and that Q and P*, representing the multiplication by x and y respectively, can 
have no more than one diagonal above the main diagonal. The converse is also 
true as will be detailed below. 

Proposition 2.1 Suppose that P and Q have the above band structure and 
that the highest diagonal of Q and the lowest of P have nonzero entries. Then 
the two following conditions are equivalent 

(i) The commutator [P, Q] is diagonal. 

(ii) There exist two polynomials of degrees d\ and d 2 respectively which we 
denote by V{(x) and V 2 ' (y) such that 

(P-y/(Q))> = 0, (Q - Kj'(P))< = , (2-25) 

where the subscripts <0 or >0 denote the lower or upper part. 

Proof. The detailed proof of this result may be found in Q. Here we just note 
that, given the band structure of the two semi-infinite matrices P and Q, the 
polynomial V{(x) may be uniquely determined from the relation 

(P - V{(Q)) ■ e = , e := (1,0,0,0,...)* (2-26) 

and its existence rests upon the assumption that j(n) ^ 0. A similar relation 
uniquely determines V^'(y). 

It may then be shown that all the relation contained in eq. (2-25) are satisfied 
by these polynomials. 

Conversely, if two polynomials V{ and V{ satisfying eq. ( [2-25Q exist, then 

[P, Q - V 2 (P)] - [P Q] = [P- Vi(Q), Q] . (2-27) 

But the LHS is upper triangular and the RHS is lower triangular (not strictly) , 
so that [P, Q] must be diagonal. Q.E.D. 



The structure (2-21), (2-22) of the two matrices P and Q means that the 
four sequences ipmip j4 > m4 > satisfy both multiplicative and derivative recursion 
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relations 



n-31 



3=0 



1 d 



y<t>n 

1 d A 



7(n - l)"0n-l + 2J ^3'( n + i)^n+J 7 
.7=0 

7(n)^ n+ i + ^2 0j{ n )<Pn-j, 
3=0 

7(n- l)0„_i + 2jaj(n + i)0„ +J - . 



(2-28) 
(2-29) 
(2-30) 
(2-31) 



From the finite recursion relations satisfied by the quasi-polynomials {ibn.(x]\ 
and {4>n(y)} follows a set of "generalized Christoffel-Darboux relations Jl^, g], 

N N 

which imply that the kernels Ku(x, x ) and Kti^V , y) may be expressed as: 



Kn(x,x ) = 



^22(2/ ,y) = 



JV-l JV 

$ (x ), a*(^; 

JV 



JV-l N 

* (w ),B*(y) 
y' - v 



(2-32) 
(2-33) 

(2-34) 
(2-35) 



where ty(x) , $(y), (y) and $ (x) are the column or row vectors of di- 

N JV 

mension (d\ + 1) and (d% + 1) defined by 

















-y(iv-i)l 


JV 


«d 2 (A0 




Q 2 (A0 


ai(A0 





A := 





« d2 (iV+l) 




Ql(JV+l) 













a d2 (iV+2) 


















a d2 (iV+d2-l) 



















-y(AT-l)l 


JV 






/3 2 (iV) 







B := 









Pi(N+l) 













j3 dl 0V+2) 


















/3 dl (N+di-l) 







JV-l 


JV- 


-1 







= [^JV-^21 ■ ■ •tV'jv]' , $(?/) = [0JY-^ 21 ■ 
JV JV 



^JVJ 



(2-36) 



JV-l 



JV-l 



£ (y) = . . . ,^ 2 _J , 1 (*) = . . ..^.(MT) 



JV JV 



The matrices A, B entering eqs. 2-35 define two pairings (which we will refer 

JV-l 

to as the Christoffel-Darboux pairings) between and $ and between $ and 



JV 



JV 



JV-l 



£ . We call these pairs dual windows. 
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The key observation is that any quasipolynomial tpj(x) can be uniquely ex- 
pressed, for any given N > in terms of linear combinations of any d 2 + 
1 consecutive basis elements ^>jv— d 2 , •-, V'iV with polynomial coefficients. We 
call this procedure folding of the space onto the window spanned by \Etjv = 
[V'jv— d 2 j ■■iV'Af]'- This is accomplished by means of the cc-recursion relations for 
the quasipolynomials in eq. (2-31), which allow us to express the (N + l)st 
quasipolynomial in terms of the d 2 + 1 preceding ones, but with coefficients that 
are polynomials in x. Iteration of this procedure defines the folding. 

In matricial form the above can be expressed as follows 



a fx) ^(x) 

N N 



<5 (x) 

N+1 



N>di . 



(2-38) 



where 



sl(x) := 

N 








-ad 2 (N) 
~/(N) 













1 

-ai(iV) (x-ao(N)) 



■y(N) 



TP) 



(2-39) 



The matrix a is invertible, since its determinant equals ad 2 {n)/^f{n) and ad 2 {n) 
can be shown not to vanish as a consequence of the relation 



(Q-V 2 '(P))< = 



(2-40) 



It will be referred to in the following as a "ladder" matrix. A completely analo- 



gous relation can also be shown for the quasipolynomials 4> n (y) (see eq. ( 2-50 ) 
below) and for the respective Fourier-Laplace transforms. 

By means of this folding, one can also express the action of any operator of 
finite band size as a matrix polynomial in x of size d2 + 1. The most relevant 
case is the folding of the operator P — — 4^ . Introducing the notation 



* := [i/iN-ch, ■ 

N 



(2-41) 



we have 



Id N 

K dx n N 



/ N N J_L N \ 

= 7 ( a (x))-^ fa +Y a (x) a (x) ■ ■ ■ a(x) (2-42) 

I N-l ' N+j-I N+j-i N I N 



where 



Pi 



dia l g[p j (N + j-d 2 ),l3 j (N + j 



da + l),.:,0j(N + j)] ,(2-43) 



7 



A r 

7 



diag [ 7 (JV - 1 - da), j(N-d 2 ),..., 7 (JV + d 2 - 1)] . (2-44) 



The corresponding statement for the <j) n 's is obtained by interchanging x and 
y, tp n and <p n , d\ and d 2 , aj and /3j etc. One obtains a similarly defined matrix 
D 2 (y) representing the action of the derivative on the quasipolynomials 0„'s. 
With the notations 



N 
N 



N 



b(y) 

N 



one finds 



'N—d\ i 



;<PN\ , 



(2-45) 



diag [ aj {N + j - di), aj{N + j - di + 1), . . . , aj (N + j)} (2-46) 
.7=0,... d 2 



7 := diag[7(JV-di-l),...,7(JV-l) 








7 (JV) 













1 

T (JV) 7 (JV) 



N > di 



(2-47) 
(2-48) 



b(y) $(y) = $ (y) 

AT AT Af+1 
Id N 

K dy n N 
:= 7 ( b rHy)+ao 

N-l 



(2-49) 



-X^aj b (y) b (y)---b(y) $. (2-50) 

J JV-W-l AT / A r 



We can repeat a similar procedure for the respective Fourier-Laplace trans- 
forms. The relevant definitions and relations are given by the following formulae 



a(ar) 



N 

h(y) 



-y(N-l) 

- ai (N+l) 
-y(N-l) 



-a d2 (N-+d 2 ) 
■y(N-l) 

y-MN) 

7(JV-1) 

-Pi(N+l) 
'((N-l) 







7 (JV-1) U 





'•• 



1 





1 







1 




e gld 2 +i[x] 



e gld 1+ i[y] 



(2-51) 



(2-52) 
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JV-1 



N N 

*a(x) , 



N-1 



N N 

*b(y) 



i d V/ \ 

i?d^ (y) 



N-1 N 

* G/)£a(tf) 



1 iV-l iV-l iV 



> di + 1 . 
N > d 2 + 1 



where 



JV TV 

Mv) ■= (b)" 



AT-l JV-1 



JV-1 JV-2 



TV 



7~+'a ~ + ^ b b • 

,N. , JV-1 ^^T 1 ^JV-lJV-2 

(a) 7 +/9o +2^ a a 



JV-j /v_i 

■ b a,- 



JV-1 
JV— 1 

0.- 



diag(a J -(JV-l),...,a J -(JV + di-l)) , 
diag(/3 i (iV-l) I ...,/3 J -(iV + d 2 -l)) 



N-jN-l 

■ a £. 



(2-53) 
(2-54) 

(2-55) 



(2-56) 

(2-57) 

(2-58) 
(2-59) 



Summarizing, we have four sequences of linear differential systems 



Size (d 2 + 1) x (d 2 + 1) 



Id N 

K dx jv jv 

1 d N-1 JV-1 JV 



Size (di + 1) x (di + 1) 



1 JV-1 JV-1 JV 

Id N 
K dy n jv 



(2-60) 



as well as the ladder relations ( 2-38 ), ( 2-49 ) , ( 2-51 ), ( 2-53 ). We have not con- 
sidered here the deformation equations, i.e. the differential equations obtained 
from infinitesimal variations of the coefficients of the potentials V± and V2 en- 
tering the measure. The complete study of these deformations is carried out in 
H . In particular it is shown there that the resulting overdetermined system of 
PDEs is compatible. Here we will only recall that the mixed system of ODEs 
and difference equations is also compatible, as implied by the following: 

Proposition 2.2 The ladder matrices a intertwine the differential systems D\ 
with different N's, i.e. 



sl(x) 

JV 



( d N 



d N + 1 

— + £>i (x) ) a(ar) 



/v 



(2-61) 



Similar statements hold for the other three sequences of ODEs and ladder rela- 
tions. 



The next proposition explains how the four sequences of systems in the Table 
are related amongst themselves by means of the Christoffel-Darboux pairings. 
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Proposition 2.3 The following relations are satisfied 



N N N N N N N N 

D 1 (x)A = ADi(x) ; D 2 (y) 1 = B D 2 (y) (2-62) 

The spectra of the two matrices Di(x) and D_\{x) (i- e -j their characteristic 
polynomials) coincide, as do the spectra of D 2 {y) and D 2 (y). 

A less apparent spectral duality also holds. Indeed it is proven in || that 

det (yl-D^x)) = cdet (xl-D 2 (y)\ , (2-63) 

where c is the ratio of the leading coefficients of the two potentials V% and V 2 . 
Notice that the two determinants involve square matrices of rank d 2 + 1 on the 
LHS and of rank d\ + 1 on the RHS. In the following section we give a simple 
derivation of a "naive" N — > oo limit of these results; namely one in which we 
treat the relevant recursion matrices as commuting. 



3 The Abelian case 

In this section we derive the spectral duality property in a particular limit N — > 
oo, K/N — 0(1). In such a limit the two matrices P and Q, while retaining their 
finite band structure, may be taken to commute because [P, Q] = — 1/K — > 0. 
In addition, we consider only the case in which the coefficients aj(n), Pj(n), 
g(n) do not depend on n: this is a stronger requirement which occurs actually 
only for certain ranges of the coupling constants. This limit is studied in the 
literature and is referred to as the "one-cut case" or the "genus 0" case j?]. ^ . 

A further simplification that is purely technical is obtained by considering 
the matrices as doubly-infinite, i.e. of size Z x Z instead of N x N. We will show 
that the statement of spectral duality in this case reduces to a classical result in 
commutative algebra, namely the computation of the resultant of two Laurent 
polynomials. 

The non-abelian case (i.e. for finite N) is detailed in || and the approach 
used there may be used to derive the result for the N — > oo case. However, we 
will present a proof here of a different nature, which can also be extended to 
the non-abelian case Q. 

The equations [P, Q] — — -hi in the limit N — > oo, K — O(N), become 
commutativity equations [P, Q] — 0. Moreover, since we are considering finite 
band matrices and we focus on the window at N, we can replace the semi-infinite 
matrices P, Q by doubly infinite matrices with the same band structure. For 
suitable scaling regimes it can be argued on physical grounds that the sequences 
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j(n), aj(n), 0k(n) actually do not depend on n (provided n = O(N)). It is 
precisely this very simple case that we want to address here. 

The pair of commuting matrices P and Q with the same band structure as 
before now just become polynomials in the shift matrix. All the matrices are 
taken to be Z x Z matrices and hence A = [J^+i] is actually invertible, the 
inverse being just the transpose A*. With this in mind we can write 



Q(A) := 7 A + a + ^ a 4 A , 7 ^0^ a d2 



di 



P(A) := 7A- 1 + 0o + M l , 7 ^ ± dl 



(3-64) 



(3-65) 



with Q and P are viewed as Laurent polynomials in A, A -1 . It is convenient to 
introduce an indeterminate A and represent Q and P as acting on the graded 
space 

Q, P : C[A, A -1 ] — > C[A, A -1 ] , (3-66) 



determined by substituting A by A in the relations (3-64,3-65). The shift matrix 



A is just multiplication by A while A* = A 1 represents multiplication by A . 
The equivalent of a window is then the linear span of di + 1 consecutive powers 
of A 

C{i) N - d2 ,...,iJ N } ~C{\ N - d2 ,...,\ N } . (3-67) 

The folding of the previous sections here reduces to a very simple expression. 
Indeed, folding the graded space W := C[[A]] onto the span of X N ~ d2 , .... X N 
simply means taking the quotient 



C[[A]] ~ C[x] <8> C[[A]] mod (x - Q(A) = 0) 



c}{\ 



N-d-2 



X N } . (3-68) 



In other words, the power X N+1 can be re-expressed in terms of the powers 
\ N ~ d2 , ...,X N using the relation x — Q(X) = 0. The equivalent of the ladder- 
matrix is just the expression of multiplication by A in the "folded" window 
C[x]{X N ~ d2 , A^}. It is defined so as to make the following diagram commu- 
tative 

C[[A]] > C[[A]] 



(p-Q(X)=0) 



C[x]{X N - d2 ,...,X N }^l 



^r-Q(A) = 0) 



(3-69) 
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In principle, a could depend on N, but it is easy to see that in fact it is repre- 
sented by the following iV-independent companion-like matrix 



a(x) — 






7 



Qj X — ag 

7 7 



(3-70) 



Similarly, we could define another folding along P by means of the following 
diagram 



C[[A]] 



A 



C[[A]] 



^/-P(A) = 0) 



(3-71) 



C{y}{\ N ~\ \ N + d ^} ^M. C[y]{X N -\ \ N + d ^} 
where b is given by 

1 ■■• 



b(y) = 






Pd 1 



1 

,§i y-Po 

7 7 



(3-72) 



In this framework the matrices D\{x) and D2{y) are simply 

di 

D x (x) :-P(a( a; ))= 7 a(x)- 1 +^/3,a(x) 

j=o 

d 2 , 

D 2 {y) ■■=Q(b(y))=ib(y)- 1 +22 a j*>(y) ■ 

3=0 



(3-73) 
(3-74) 



The previous statement about spectral duality now translates into the iden- 



tity 



dct(?yl - Dx{x)) oc det(xl - D 2 {y)) 



(3-75) 



We will show that both determinants are in fact the resultants (w.r.t. A) of the 
two Laurent polynomials Q(A) — x and -P(A) — y. The proof is actually quite 
standard for polynomials and here we just adapt it to the situation with Laurent 
polynomials (see e.g. JTo|). This amounts to studying the following embeddings 



C{A 



N-d 2 



,x N } c{\ N - d *-\ 



,X N+dl }, 



N-i ^w+di-ij £{X N ~ d2 ~ 1 X N+dl } 



(3-76) 
(3-77) 
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where x and y are treated as parameters of the embeddings. Let us denote by 
W, W_, U the three vector spaces 

W := C{X N - d2 ,...,X N } , W:=C{X lf - 1 ,...,X N+dl - 1 }; 
U := C{X N - d2 -\...,X N+dl } . (3-78) 

The above embedding may be combined into a single map 

W ®W_ — — — — — — — — lu (3-79) 

The two parts of this map give spaces generically transverse as x and y vary If 
they are not transverse for a given pair (x,y) 7 this means that 

3w e W, 3w £W_ such that i»/0^i», (P-y)w = (Q-x)w € U . (3-80) 

Taking the quotient of this relation by the relation Q(X) — x = or P(X) — y = 
gives rise to the relation 

(D 1 (x)-y)w = (3-81) 
(D 2 (y)-x)w = , (3-82) 

which means that y is an eigenvalue of D\(x) and x an eigenvalue of D-i[y). 



Conversely if either of the two equations (3-81) (3-82) holds, say the first, for 



nonzero vector w, this means that there exists w € W_ such that 

(P - y)w = (Q- x)w . (3-83) 

Notice that (Q — x)w cannot be zero since the map P — y : W — > U is injective 
for all y and so (P — y)w ^ 0. The same result follows if we start from eq. 
( ||f ). This proves that the embedding is not transverse if and only if x is an 



eigenvalue of D2(y) which is equivalent to y being an eigenvalue of D\(x). 

The condition of transversality amounts to the nonvanishing of the deter- 
minant of the embedding (in any fixed basis). It is easy to see that such an 
embedding is represented by the Sylvester matrix 
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of the two Laurent polynomials, whose determinant A(x, y) equals the resultant. 
A simple counting of degrees and inspection of the highest powers in x or y shows 
that 

a d2 j dl det(yl - D 1 (x)) = A(x,y) = f3 dll d2 det(xl - D 2 (y)) . (3-85) 

which defines the spectral curves as the non-transversality locus of the embed- 
dings. The intersection of the two embeddings on this spectral curve is (gener- 
ically) one-dimensional and projects to the eigenvectors of D\(x) and Z?2(y)- 

While this is very simple, and just a reformulation of standard algebraic 
results in this abelian setting, a very similar approach can also be used to prove 
spectral duality for the pair D\ N (x) and D 2 (y) in the finite N setting, in 
which the matrices P and Q do not commute. A refinement and elaboration on 
this theme also leads to the other results of Q in a more elegant and compact 
form Q , such as the compatibility of the deformation equations in the coupling 
constants of the potentials V\, V2 which, in particular imply the invariance of 

N N 

the generalized monodromy of the operators d x + Di(x) and d y + Dzty). This 
defines a sort of "noncommutative resultant" for finite band matrices whose 
properties will be developed in a subsequent publication. 
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